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Hysteretic I-V (current-voltage) is studied in narrow Al nanowires. The nanowires have a cross 
section as small as 50 nm 2 . We focus on the retapping current in a down-sweep of the current, at 
which a nanowire re-enters the superconducting state from a normal state. The retrapping current 
is found to be significantly smaller than the switching current at which the nanowire switches into 
the normal state from a superconducting state during a current up-sweep. For wires of different 
lengths, we analyze the heat removal due to various processes, including electronic and phonon 
processes. For a short wires 1.5/xm in length, electronic thermal conduction is effective; for longer 
wires 10/xm in length, phonon conduction becomes important. We demonstrate that the measured 
retrapping current as a function of temperature can be quantitatively accounted for by the self- 
heating occurring in the normal portions of the nanowires to better than 20 % accuracy. For the 
phonon processes, the extracted thermal conduction parameters support the notion of a reduced 
phase-space below 3-dimensions, consistent with the phonon thermal wavelength having exceeded 
the lateral dimensions at temperatures below ~ 1.3K. Nevertheless, surprisingly the best fit was 
achieved with a functional form corresponding to 3-dimensional phonons, albeit requiring parameters 
far exceeding known values in the literature. 



Understanding the dynamics of ultra narrow supecon- 
ducting (SC) nanowire wires is an active area of in- 
vestigation [1H15|. A significant focus is the so- called 
1-dimensional (ID) limit, delineated by the condition, 
(w, h) < £, where w is the width and h the height of the 
nanowire, and £ the superconducting coherence length. 
Investigations of the behavior under current-biasing not 
only elucidate the conditions and limitations for the cur- 
rent carrying capabilities, as well as the process of recov- 
ery back into the superconducting (SC) state after driven 
normal by an excessive current, but also potentially lay 
the foundation and pave the way for the development of 
novel devices, such as a current- Josephson effect devices 
0, or qubits 0. 

In this work, we report on measurements carried out in 
ultra narrow Al nanowires with a cross section as small 
as 50 nm 2 . The three nanowires studied have a widths 
and heights ranging between 7-10 nm, and lengths of 1.5 
//m (wire SI) or 10 jum (wires S2, S4). These nanowires 
are exceedingly uniform in their cross section, as indi- 
cated by their ability to carry sizable current before being 
driven normal, where the current density is nearly iden- 
tical to co-evaporated 2d films. In a previous work, the 
behavior of the switching current I s during an up-sweep 
of the current was investigated Q. There, it was found 
that heat deposited by phase-slips- transient temporal- 
spatial events during which the superconducting phase 
fluctuates and changes by 27r over a distance of order 
£, while the core region goes normal-leads to a thermal 
runaway, driving the entire nanowire into a normal state 
from the SC state. 

Here we focus on the down-sweep retrapping current. 
The retrapping current I r is found to be significantly 
smaller than the up-sweep switching current I s , and can 
be a much as a factor of 20 smaller. The history de- 



pendent current-voltage (IV) relation exemplified by the 
disparate behaviors in the up- and down-sweep is ubiq- 
uitous, despite the fact that based on the criteria nor- 
mally applied to SNS (superconductor-normal metal- su- 
perconductor) bridges, the nanowires should be in the 
heavily over-damped regime in its dynamics 0, 0, [H| ■ In 
MoGe nanowires of widths ~ lOnm, Tinkham et a/.[(| 
performed a heat flow analysis, and ascribed the retrap- 
ping behavior to self-heating. Our work bears similarity 
to that work, but our SC nanowires are in a different 
regime, where fcpZ ~ 60 3> 1, rather than being close to 
1 in their case. Here fc^ is the Fermi wave-number, and 
I is the mean-free-path. Moreover, their nanowires were 
suspended freely, while ours are deposited onto a nar- 
row, 8 nm-wide InP ridge (Fig. [lja)), and are thus in 
thermal contact with an underlying substrate. Further- 
more, our analysis differs from theirs in the form of the 
heat flow equations. Based on our analysis, we rule out 
under-damping as the cause of the hysteresis, in agree- 
ment with recent results in submicron SNS bridges [la ]. 

To lay the framework for understanding the behav- 
ior of nanowires, the Josephson junction can serve as a 
starting point. There, the free energy landscape under 
current bias is described by the tilted washboard poten- 
tial, shown in Fig. Ufa) [ljj]. This same scenario is 
also applicable to ID SC nanowires [ill [l2|, Josephson 
junctions are classified within a Resistively and Capaci- 
tively Shunted Junction (RCS J) model as either under- or 
over-damped, depending on whether the quality factor, 
Q = \/2eI c C /hR, is greater or less than 1. Here, I c is the 
critical Josephson current, C is the junction capacitance, 
and R the junction normal state resistance. When under- 
damped Josephson junction is driven over the free-energy 
barrier out of its meta-stable minimum, the SC phase 
keeps running downhill as there is insufficient damping 
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FIG. 1: (a) Schematic of the Al superconducting nanowire 
device on a narrow InP ridge template. The ends of the 
nanowire are connected to large, electrical measurement pads. 
The pads can either be in the superconducting state, or driven 
normal by a magnetic field, (b) Top view of the nanowire ge- 
ometry, and the layout used in the heat flow model discussed 
in the text. 



to retrap the phase in a lower energy local minimum. 
A consequence is that a hysteretic current-voltage (IV) 
relation, where the up-sweep and down-sweep branches 
do not overlap. In contrast, in an over-damped junction, 
the phase moves diffusively between adjacent minima, 
and hysteresis is often not present |2ll423j . 

The estimated Q for our nanowiresis in the range of 
~ 0.01, far below unity, and the nanowires are ostensi- 
bly in the severely over-damped limit. This estimate is 
relevant when the nanowire device is in the S-NW-S con- 
figuration, where S refers to each of the two large metal- 
lic measurement pads when in the SC state, and NW 
denotes the nanowire. It also provides a reasonable es- 
timate in the N-NW-N configuration, when the pads are 
driven normal, but the ambient temperature is below the 
nanowire SC transition temperature T c . In this case, the 
nanowire itself breaks up into alternating SC and normal 
segments, whether during an upsweep or a downsweep of 
the current. In the former case, the nanowire is overall 
in the SC state, but during a phase-slip, the phase-slip 
normal core acts as the normal region. In the latter, the 
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FIG. 2: (a) Tilted-washboard free energy landscape for a 
Josephson junction under current bias. A phase slip occur- 
ring between adjacent minima is indicated. A similar scenario 
occurs in a ID superconducting nanowire. (b) Hysteretic IV 
curves for nanowire S2 at several temperatures. Arrows in- 
dicate direction of current sweeps. The upsweep switching 
current I s is significantly larger than the down-sweep retrap- 
ping current I r . 



central portion of the nanowire is normal due to heating, 
while the regions closer to the pads are in the SC state 
(Fig. 1(b)). 

Nevertheless, despite the over-damping, hysteretic IV 
curves are ubiquitous, as can be seen for wire S2 in Fig. 
[Ha). In fact, the ratio of I r to I s can be as small as 
~ 1/20. For example, in nanowire S2 at T ~ 0.3A", I r ~ 
0.19/zA, while I s ~ 4/iA These observations motivated 
us to investigate the retrapping current systematically, 
as a function of the temperature and wire length, and to 
perform a detailed heat analysis to establish self-heating 
as a cause of the substantially reduced I r below the value 
of the upsweep / s . 

Our devices were fabricated using a template method. 
The template is a narrow, 8nm wide InP ridge, formed 
by differential etching on the cleaved (110) crystallo- 
graphic plane of a molecular-beam-epitaxy (MBE) grown 
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InGaAs-InP crystal, where the growth direction is (001). 
The geometry of our devices is depicted in Fig. QJa). 
The details of the fabrication procedure is described in 
a previous work [20j ]. The nanowire resides on the nar- 
row InP ridge and is thus thermally connected to the 
large semiconductor substrate through the narrow ridge. 
The nanowire is electrically connected to large metal- 
lic measurement pads on its ends. Therefore, for heat 
removal, thermal conduction both in the lateral direc- 
tion along the nanowire, and vertically through the InP 
ridge via phonon processes must be considered. The IV 
measurements were carried out in a extremely carefully 
shielded apparatus to minimize unwanted environmen- 
tal interference, such as external noise (e.g. from nearby 
radio stations) conducted down the electrical cables, or 
Johnson-Nyquist noise from resistors within the electri- 
cal measurement circuitry. In particular, Thermax ca- 
bles with the ability to remove high frequency noise is 
employed where possible, as well as low-temperature RF 
filters. The devices are also enclosed in metal cans with 
all openings plugged with conductive tape or metal mesh. 

In Fig. EJa)-(c), we present the measured I r for 
the three nanowires SI (lOnm x lOnm x 1.5/im), S2 
(9.3nm x 9.3nm x lQ/im), and S4 (7nm x 7nm x 10^m), 
respectively, as a function of temperature. Immediately 
apparent is the substantial difference in the magnitude 
of the retrapping current I r for the short SI, when com- 
pared to the longer S2 and S4. For SI, I r = 0.99fiA at 
T = 0.3K, while it is 0.19^,4 and 0.117/M, respectively 
for S2 and S4. Thus, the value of I r in SI is roughly 5 
- 8 times that in S2 or S4. The results for S2 and S4 do 
not depend on whether measurements were carried out 
in the S-NW-S or the N-NW-N configuration. On the 
other hand, for SI, due to the larger current level and as- 
sociated excessive heating, it was necessary to drive the 
measurement pads into a normal state, into the N-NW- 
N configuration. By driving these pads into the normal 
state using a small magnetic field B — 0.1T, the now 
normal pads can act as good thermal anchors, keeping 
the temperatures of the ends at the ambient tempera- 
ture T a . For the longer S2 and S4, the smaller current 
level means that a much lower amount of heat needs to 
be carried out through wire ends; thus improved thermal 
anchoring was not needed. For SI, on the other hand, 
keeping the pads SC (S-NW-S configuration) reduced I r 
to ~ 0.55 [i A from ~ lfiA as heat removal becomes more 
difficult due to poor thermal conduction capability of the 
SC pads. At the same time, instabilities arise in the tem- 
perature profile along the nanowire, leading to very noisy 
data below 0.8K with I r fluctuating as much as 0.1 [iA 
between adjacent data points. 

The configuration with normal electrical pads is rele- 
vant for our data in Fig. 3. In the analysis which fol- 
lows, we will focus on this configuration. The retrapping 
process returns a nanowire into the SC state during a 
down-sweep of the current, /, at an ambient lattice tem- 



perature, T Ql below the zero current critical temperature, 
T C (I — 0). At large /, most of the nanowire remains nor- 
mal due to self-heating, which raises the local temper- 
ature above the switching temperature at that current, 

During a current downsweep, the temperature is posi- 
tion dependent along the nanowire. In the N-NW-N con- 
figuration, each of its ends is connected to a large, normal 
metal electrical measurement pads, anchored at T a . Be- 
cause T a < T C (I = 0), and is in fact below T S (I = I r ), 
the end regions are in the SC state, aside from a short 
proximity region, l prox ~ lOOnm in length, immediately 
adjacent to each normal pad (Figs. Q~tb)). Here, the 
T s (I) is the upsweep switching temperature at current J, 
On the other hand, as long as / exceeded I r , the center 
of the nanowire is above T s (I) and is thus in the normal 
state. When / is reduced down to I r , a blockage prevent- 
ing the central region to be cooled is suddenly removed, 
allowing the cooling to propagate all the way to the cen- 
ter. This blockage removal is described in detail in what 
follows. 

The upsweep switching temperature for a given current 
T S (J) expresses the same relation as the switch current 
as a function of temperature I S (T), but viewed in re- 
verse. For our nanowires, these were reported in Ref. 8 
0. It is worthwhile to point out that the value of I s at 
a given temperature is not unique, but depends on the 
upsweep ramp rate of the current. Conversely, T S (I) is 
also dependent on the ramp rate. This is due to the fact 
that the switching is caused by phase-slip events 043, 
and thus the relative rates of the phase-slip generation 
to current upsweep matters. The slower the upsweep 
current ramp rate, the longer waiting time is available 
for phase-slip events to take place within a given interval 
in current. This increases the probability for switching 
within that interval and reduces the magnitude of the 
current at which switching takes place. But because the 
rate of phase-slip generation is exponentially dependent 
on the current, the dependence of I s on ramp rate is 
weak, and is approximately logarithmic. 

To estimate I r , it is necessary to determine the position 
dependent temperature, T(x). Both the electronic and 
phonon thermal conduction mechanisms need to be con- 
sidered. The temperature range of interest is low com- 
pared to the lattice DeBeye temperature Qd ~ 300K 
and the temperature dependence of the phonon ther- 
mal conductivity takes a power-law form, reflecting the 
phonon density of states. At these low temperatures, 
0.2K < T < 1.3K, the value of the phonon conductiv- 
ity is considerably smaller than the electronic thermal 
conductivity. The only exception is in the SC regions 
when T(x) < 0.35^. Thus, in the absence of a very 
large temperature rise phonons can only carry away a 
relatively small amount of heat. Whereas electronic con- 
duction requires the heat to exit the ends of the nanowire, 
the phonon conduction goes through the short InP tern- 
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3: Solid squares-Data for the retrapping current 
a function of temperature for nanowires (a) SI 
(b) S2 (9.3nm x 9.3nm x 10/im), 



FIG 

/ r as 

(lOnm x lOnm x 1.5/xm) 
and (c) S4 (7nm x 7nm x 10/im). For SI in (a), the solid 
curve is a fit based on the heat flow model describe in the 
main text, where only the electronic thermal conduction is 
considered. The Wiedermann-Franz constant deduced from 
the fitting is L a = 3.65 ± 0.15 x 10- 8 WQ/K 2 . Vertical bars 
indicate certainty of the fitting curve, due to a combination 
of the uncertain in the parameters, including uncertain in the 
superconducting transition temperature T C (I) of ±2%. For 
S2 and S4 in it is necessary to include both electronic and 
phonon conduction. Various rate limiting phonon processes 
were considered, including electron phonon relaxation and 
Kapitza boundary resistance phonon conduction (see text). 
The uncertain of the phonon parameters is ±10%. 



plate ridge (in height), on which the nanowire resides. 
The nanowire length is microns /zms while the ridge is 
only 30nm in height. If the wire is long, the electronic 
mechanism will become much less effective, and phonon 
conduction must be included as well. 

For the short wire SI, electronic conduction over- 
whelmingly dominates. For the long wire S2, both elec- 
tronic and phonon thermal conduction must take place 
side by side. Phonon conduction through the InP ridge 
takes place via several steps: (a) electron-phonon energy 
relaxation within the aluminum nanowire, (b) conduction 
through the aluminum- InP boundary, and (c) conduction 
through the 8 nm wide, 30 nm tall InP ridge. Below the 
ridge, the energy is dissipated in the very highly con- 
ductive GaAs bulk material. Thus the base of the InP 
ridge can safely be assumed to be held at the ambient 
temperature, T Q . 

One additional mechanism of heat removal takes place 
through the liquid He 3 , which surrounded the nanowires 
in the set of measurement on SI and S2. However, this 
channel appears less important. Additional data for wire 
S4 obtained in the dilution refrigerator, in which the sam- 
ples are in vacuum and thus there is no liquid surrounding 
the nanowire, yielded a retrapping current which can be 
accounted for in a similar manner as S2, using electronic 
and InP ridge thermal conduction only. Note both wires 
S2 and S4 are 10/im in length. 

We divide a nanowire into two symmetric halves of 
length L/2 each, where L is the total wire length, and 
consider the right half, where < x < L/2, and T{x = 
L/2) = T a , as shown in Fig. [ljb). When slowly down- 
sweeping the current /, we assume that the nanowire is 
in the SC state at a position x, if T{x) < T S (I), but is in 
a normal state if T(x) > T S (I). 

We begin by considering the short wire SI (1.5/zm in 
length) and only include electronic conduction. Phonon 
conduction alone will remove ~ 3% of the heat generated 
by heating, and will be neglected. The diffusion equa- 
tions must account for three regions: (a) the central re- 
gion for < x < Xb, which is normal for I > I r and has a 
resistance per unit length of Rn / L ~ 0.33 — 0.82fcf2//im 2 
[H and thus self-heats; (b) the SC segment for x\, < x < 
L/2 — Iprox, which nearly does not self heat, but must 
conduct the heat generated by the central normal seg- 
ment, and (c) the short proximity region adjacent to the 
normal metal pads for L/2 — l prox < x < L/2 which is 
approximated as a normal region. More precisely, the SC 
region close to the SC-normal boundary Xb heats slightly 
due to occasional phase slips as its temperature is just 
below T S (I) > 1.1K; the proximity region generates self- 
heating, and in addition must conduct through it the 
heat of the central normal region as well. For the normal 
metal regions, the steady state heat diffusion equation is 
given by: 
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FIG. 4: Plot of the thermal heat capacity ratio, use I 'kn, be- 
tween the normal and superconducting states, as a function 
of reduced temperature, T/T c . The line indicates the approx- 
imation used in the calculations. 



where kn = L LT/(RNA nw ) is the electronic 
Wiedemann-Franz electronic thermal conductivity, L Q is 
Lorenz number determined from fitting, and A nw = wxh 
the nanowire cross sectional area. In the SC region, it is 
replaced by: 



« --r{KscA nw — ) 
ax ax 



(2) 



The equations in the three regions are supplemented by 
boundary conditions at the junctions. The junction be- 
tween the center normal region and the SC region takes 
place at x — xt, and between the SC and short prox- 
imity regions, at x = L/2 — l pr0 x- The temperature is 
continuous across each junction, and the heat flow is 
identical immediately to the left and right. Lastly, we 
have T(x = L/2) = T D . Note that the forms of these 
equations differ from those used in Ref. 6. There, the 
variation of the thermal conductivity kn or ksc with 
position, through their dependence on temperature T, 
was not accounted for Q. 

At a given temperature below T C (I), the electronic 
thermal conductivity of the SC segment at a position x, 
with a temperature T{x), is related to the corresponding 
normal metal W-F thermal conductivity at that temper- 
ature, by the ratio r = ksc/kat:[3] 



71-2 A.76Tc(/)/r(x) cosh?(x/2) 



dx. 



(3) 



This integral evaluated numerically is shown in Fig. 
IH For the purpose of determining the position depen- 
dent temperature at a given /, this ratio is approxi- 
mated by a linear form as indicated in the figure for 



0.2 < T{x)/T c {I) < 1: 

r « lA25[T(x)/T c (I) -0.2] 



(4) 



We next describe how the cooling blockage can be re- 
moved, and determine the condition for this to occur. To 
do so, we consider the SC region and fix the current at /. 
The SC-normal boundary occurs at x&, which is deter- 
mined by equating the temperature at Xb, T(xb), to the 
switching temperature at that current T s (J) . Momentar- 
ily treating Xb as a variable, T[xb) attains its maximum 
value at a critical value Xb = x c ; in the simplest approxi- 
mation x c is independent of / as will be seen below. For a 
large current /, the normal region is large and the actual 
Xb exceeds x c (Fig. EJ. As / is reduced down to I r , the 
normal region shrinks and Xb becomes equal to x c . Here, 
T{xb = x c ) takes on the maximum possible value at I r , 
since xj, — x c , and is equal to T s (I r ). Further reducing 
/ to just below I r , T S (I < I r ) will slightly increase from 
T s (I r ), while at every x, T(x) will slightly decrease due 
to reduced heating. The decreased maximum tempera- 
ture at x c , T(x = x c ), can no longer reach the increased 
T S (I < I r ). The boundary will become unstable, and 
will propagate toward the center at x = 0. Starting from 
the inital boundary at Xb = x c , more and more of the 
normal region will fall below T S (I) and become super- 
conducting, as the shrinking normal region generates less 
and less heat, until the entire wire is cooled. The three 
cases are depicted in Fig. [5j 

The value for x c can be deduced from the steady-state 
heat diffusion equation in the SC state, Eq. 2. Neglecting 
the short proximity region adjacent to the pad, x c — L/A 
(0.5(L/2)); accounting for the proximity region of length 
Iprox modifies this to x c ~ {L/2 — lp rox )/2. For illustra- 
tive purposes, let us determine x c in the absence of the 
proximity region. Focusing on the SC region at its border 
with the normal segment, x — Xb, twice integrating the 
diffusion equation and matching the boundary conditions 
yields for the left- hand- side (LHS): 



LHS = (I 2 R N -)(L/2 - x), 



(5) 



which is maximal for x — L/A for fixed /. At x = Xb 
the factor I 2 Rjy^ represents the heating power gener- 
ated by the normal region where < x < Xb- Equating 
the LHS to a twice-integrated right-hand-side and solv- 
ing for T{xb) thus yields the highest temperature at the 
normal-SC border when Xb = x c = L/A, where the LHS 
is maximal. 

Using the approximate form of the thermal conductiv- 
ity ratio r between the SC and normal states given by 
Eq. 4, the diffusion equation can readily be solved an- 
alytically. The solution yielded the position dependent 
temperature profile shown in Fig. |6ja). The kink at 
the SC-proximity boundary is an artifact of our model, 
where the proximity region is approximated as a nor- 
mal region. A more accurate model would require solv- 
ing the Usadel equation, which is expected to yield a 
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FIG. 5: The temperature profile versus position x at ambi- 
ent temperature T , with a corresponding retrap current I r . 
Three different current levels are depicted: I > I T , I = I T , 
and I < I r . The critical value of the SC-Normal boundary, is 
given approximately by x c = 0.5(L/2 — l pl0 x)- The proximity 
region lies to the right of x = (L/2 — l prox ), ending at the 
electrical pad at L/2. The SC segment lies immediately to 
the left. For / > I r , the position of the SC-Normal boundary, 
Xb, is given by the intersection of T(x) and the corresponding 
switching current T„(I). The horizontal dotted line depicts 
T s (I r ), while the short light solid lines correspond to I < I r 
and I > I r with [T S (I > I r )] < [T S (J P )] < [T S (I < I r )]. For 
I > I r , xt, > x c . For / = I r , xt = x c . For I < I r , indicated 
by the light solid curve joined onto the light dotted curve, the 
dotted "normal" region is unstable; the region immediately 
to the left of x c actually falls below T a (I < I r ) and will go 
superconducting, leading to a propagation of the SC-Normal 
boundary toward the center at x — 0. 



rounding and smoothing of the kink. See Ref. 24 for 
an example of the rounded voltage profile in the proxim- 
ity region. The fit to the I r as a function of tempera- 
ture T is presented in Fig. [3Ja) , using a Lorenz number 
L = 3.65 x 1Q- 8 WQ/K 2 , somewhat higher that the the- 
oretical value of 2.45 x 10~ 8 W£1/K 2 . Viewed in another 
way, forcing L Q to take the theoretical value, our model 
would predict a low temperature I r ~ 0.8/iA, rather than 
the 1/j.A we observed. Reconciling this discrepancy may 
require the development of more sophisticated analysis 
using the Usadel equation, while incorporating heating 
and a position dependent temperature at the same time. 
Despite the discrepancies, the overall behavior and mag- 
nitude (within 20 % accuracy) are captured in our sim- 
plified model. 

For the 10/j,m wires S2 and S4, electronic thermal con- 
duction alone is not sufficient to support the measured I r , 
despite its smaller value (by a factor of 5 - 8 at 0.3K). It 
is necessary to include phonon conduction. We assume a 
power-law temperature dependence in the heat exchange 
rate. At every position, the heat-removal linear power 
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FIG. 6: Temperature profile versus position x for different 
ambient temperatures. The temperature is stepped by 0.1 K 
between successive curves-(a) nanowire SI, (b) nanowire S2, 
and (c) nanowire S4. The kink at x = 0.79(L/2) in (a) is an 
artifact of modeling the proximity region as a normal metal, 
with a sharp boundary with the SC region to its left. 



density is set proportional to T 1 — TJ , where 7 may be 
varied to reflect the limiting process in the phonon con- 
duction discussed above. The steady-state heat diffusion 
equation takes the forms: 
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and 

-KTr-Wv-^so^), (7) 

for the normal and SC segments, respectively. For the 
longer wires S2 and S4, the smaller I r values allow the 
short proximity region adjacent to the normal-metal pads 
to be neglected. 

This form for the heat-removal power-density-per- 
unit-length assumes that phonon conduction along the 
wire direction is weak-a reasonable assumption since 
the thermal conductivity is small, and since the wire 
is much longer than the height of the InP ridge (i.e. 
10 /jm 3> 30 nm). It also assumes that the environment 
to which the energy is dissipated is well-anchored at 
the ambient temperature. For instance, if the electron- 
phonon relaxation is the limiting step, it is reasonable to 
assume that the combination of Kapitza and InP ridge 
phonon conduction is sufficiently large that the phonon 
temperature within the nanowire is maintained at T Q . On 
the other hand, if Kapitza boundary resistance limits the 
rate of heat removal via the phonons, then it is reason- 
able to assume the phonons within the InP ridge are at 
T . 

Several choices are possible for 7. For electron-phonon 
energy relaxation via 3-dimensional (3D) phonons, 7 = 5, 
while for Kapitza boundary resistance, 7 = 4. However, 
the small lateral dimensions (width and/or height) in ei- 
ther the nanowire or InP ridge should render the phonons 
reduced below 3D, down to ID and 2D, respectively. This 
occurs because the thermal phonon wavelength A p /j/2, 
exceeds the lateral dimension(s) for T < 1.3K (more 
precisely, A p /j/2 > (w,h)). Note that the width of the 
InP ridge is wi n p = 8nm while the width/height of the 
nanowire w or h ~ lOnm. By examining the fitted coeffi- 
cient in front of the power-law term, one seeks to exclude 
various possibilities. This coefficient (3 naturally depends 
on the limiting mechanism. For example, in the case 
of electron-3D phonon relaxation, /3 = 'E e -3D- p f l A nw , 
where T, e -3D-ph is the 3D energy relaxation rate con- 
stant and A nw the nanowire cross sectional area. For 
Kaptiza boundary resistance limited thermal conduc- 
tion with 3D InP phonons, /3 = <JK,3D~phWinP, where 
&K,3D-phT is the Kapitza boundary conductance at 
temperature T. 

These highly nonlinear equations were solved approx- 
imately by numerical methods, yielding the position 
dependent temperature profiles shown in Figs. IU(b) 
and (c), respectively for S2 and S4. The fitting to 
versus temperature yielded the curves in Figs. Eljb) 
and (c), with L D fixed at the value from wire SI, of 
3.65 x 10~ 8 WQ/K 2 . The critical value x c was found 
to shiftly slightly toward the center, to 0.23L (0.46(L/2) 
rather than L/4 (0.5(L/2)). The best fit is for 7 = 5, cor- 
responding to the electron-3D phonon energy relaxation 
as the limiting step. Fits of slightly lower quality can 



be achieved for 7 = 4 or 3. On the other hand, the nu- 
merical values for the coefficient j3 yielded values for the 
parameters, which point to Kapitza boundary resistance 
at the Al nanowire-InP interface as the limiting path to 
phonon conduction. 

For 7 = 5 the extracted electron-3D-phonon relaxation 
rate 'SeSD-ph for S2 is ~ 6x as large as the established 
value - 2 x 10 9 W/m 3 K 5 \v^ and is - 11 x for S4. The 
values are thus inconsistent. We are forced to consider 
the possibility that the Al phonons are reduced in dimen- 
sions down to ID. An enhancement of ~ 6_o/T(a/w) ~ 9 
can be expected per dimension reduced, where a is the 
lattice constant, yielding a factor ~ 81, far larger than 
the measured enhancement! Instead, within this sce- 
nario, one expects the limiting step to be the Kapitza 
boundary resistance between the Al nanowire and the 
InP ridge. 

For Kapitza boundary limited conduction, one may ex- 
pect 7 = 3 rather than 7 = 4 due to reduced-dimension 
2D InP phonons, despite the poorer quality fit. An en- 
hancement factor of ~ 9 over the known 3D value should 
be present from the reduction in dimension by 1 . As a ref- 
erence, we use (JK,3D-ph ~ 

2QW/m 2 K i 1 obtained for Au 
on GaAs rather than Al on InP [25| . Forcing 7 = 4 yield 
the parameter value crK,3D- P h ~ l90W/m 2 K 4 , which is 
10 times the reference value. Instead, after conversion of 
the reference value to account for 2D phonons, the en- 
hanced value of <JK,2D-ph ~ 280W/m 2 K 3 is consistent 
with the fitted values of 300W/m 2 K 3 and SZOW/m 2 ^, 
respectively, for S2 and S4 (with 7 = 3). To make certain 
this picture of Kapitza boundary limited conduction is 
consistent, we need to ensure the ridge phonon thermal 
conduction is larger. An estimate of the phonon ther- 
mal conductivity of the InP ridge itself yields a lower 
bound of 700W/m 2 K 3 , corresponding to the case of a 
very short, ridge-width limited phonon mean free path 
~ 8nm. This rules out phonon conduction in the InP 
ridge as the limiting step, as required. Finally, as a refer- 
ence, we estimate the conductivity in the absence of any 
electronic contribution. This yield a value roughly double 
the above values for the 2D phonon Kaptiza boundary 
resistance coefficient. 

It is worthwhile to re-emphasize the evidence for re- 
duced dimension Kaptiza boundary conduction as the 
limiting step, based on a direct comparison of the fitting 
parameter values for S2 and S4. Whereas, Y, e -3D-ph for 
electron-phonon relaxation limited heat removal shows a 
discrepancy between S2 and S4 of a factor 1.2/2.2 ~ 1.8, 
the values for the reduced dimension o~K,iD-ph for 2D InP 
ridge phonon, Kapitza boundary limited thermal conduc- 
tion are within 10% of each other! This, in conjunction 
with the discrepancy with the known reference value for 
^•e-3D-ph, helps establish the Kapitza boundary resis- 
tance limited scenario. 

Based on the detailed analysis presented in this work, 
we establish that in the retrapping process, the longer 
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wires S2 and S4 require phonons to contribute to heat 
removal, in addition to the electronic thermal conduc- 
tion, while for the short SI, electronic conduction alone 
is sufficient. The reasonable fits using sensible parame- 
ters demonstrate that it is possible to achieve an under- 
standing of the heating-induced hysteresis for nanowircs 
SI, S2 and S4, based on heating within the normal re- 
gions, while at the same time account for the observed 
differences. 
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